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We analyze high-field current fluctuations in degenerate conductors by mapping the electronic 
Fermi-liquid correlations at equilibrium to their semi-classical non-equilibrium form. Our resulting 
Boltzmann description is applicable to difi'usive mesoscopic wires. We derive a non-equilibrium 
connection between thermal fluctuations of the current and resistive dissipation. In the weak-field 
limit this is the canonical fluctuation-dissipation theorem. Away from equilibrium, the connection 
enables explicit calculation of the excess "hot-electron" contribution to the thermal spectrum. We 
show that excess thermal noise is strongly inhibited by Pauli exclusion. This behaviour is generic 
to the semi-classical metallic regime. 
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I. INTRODUCTION 



. In this paper we address a technologically important open problem: non-equilibrium noise in strongly driven 
^ ■ degenerate conductors. Nano-fabrication has made possible a variety of refined measurements of transport and 
I ' noise, for many different structures at sub-micron dimensions Alongside the experiments there has been much 

TiS theoretical activity . 
^ A thriving topic is the behaviour of current fluctuations in diffusive vifires. Typically, this concerns structures shorter 

than the bulk inelastic mean free path but still much longer than that for elastic scattering. They are in an operating 
region where randomness of the carrier motion prevails. This is our regime of interest. 

Two theories have come to the fore as methods of choice for describing mesoscopic transport. These are inherently 
weak-field models, predicated upon exclusively linear forms of transport analysis. One technique (Landauer-Biittiker) 
is based on coherent quantum transmission This has been adapted to fluctuations and noise by Khlus [p^ , 

Lesovik Beenakker and Biittikcr p^ , p^ , Martin and Landauer and many others [p|, r 7| . Another approach 
(Boltzmann-Langevin) uses stochastic transport equations H, reduced to a diffusive model|l8| . The same phe- 



T— I ■ nomenology has since been applied to fluctuations and noise by Nagaev |15| and de Jong and Beenakker [p|Jlq]. 

■ Although these mesoscopic-noise methodologies are markedly distinct, both agree on their paradigm. They view a 
mesoscopic wire as a random assembly of individual elastic scatterers, in a bath of free carriers whose propagation, 
impeded by the scatterers, must be regarded as strictly diffusive and must be calculated as such [ p^ . For an exhaustive 

, survey of diffusive noise theories we cite the recent review of Blanter and Biittiker [pT| as well as the earlier one of de 

a Jong and Beenakker j^] . 
The leading diffusive methods for noise also share a number of difhculties. Not least among these is the issue of 
. conformity with the fluctuation structure of charged Fermi liquids p9| . This and other basic problems are analyzed 
^ ■ in detail in references |20| and 

O i Clean, truly first-principles solutions certainly exist for non-equilibrium noise. Monte Carlo simulations are very 
O ■ well established, as witness those of the Lecce group , even if high-field Monte Carlo is still rare for metals |^,|^ . 
^ , In the non-degenerate case an analytical, self-contained and computable theory of noise has been formulated by 
JJ"! ■ Korman and Mayergoyz Their approach is strictly kinetic and free of superfluous phenomenological props. In 
rN ' philosophy it is akin to Green- function models for fluctuations, such as Stanton and Wilkins' [^ , ^ . 

' In developing a microscopically consistent account of noise, there are cogent reasons to stay within the traditional 
" " kinetic canons rather than embrace novel diffusive doctrines. The chief reason, perhaps, is sheer technological need; 
device designers can scarcely afford to be hobbled by phenomenologies whose congenital linearity denies any access 
to the vital high-field region. The small scale of modern device structures means that they are routinely driven into 
non-linear response . 

To illustrate this point we estimate the range of validity for linear diffusion. In the weak-field limit the Einstein 
relation p^ , |20| ], or drift-diffusion equivalence, underpins diffusive transport. Roughly speaking, drift-diffusion equiv- 
alence breaks down (and linear diffusive transport with it) when the energy gained in drift mediated by inelastic 
scattering exceeds the energy scale for diffusion mediated by elastic processes. In a short metallic wire this means 
that (eF/L) min{L, Lin} ^ ^vp/Ld, where V is the driving voltage, L is the sample length and Lin is the (bulk) 
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inelastic mean free path at the Fermi surface. Similarly L^i is the elastic mean free path (wp is the Fermi velocity). 
Note that for Li„ > L the effective inelastic path becomes the sample length, since dissipation in the bounding leads 
is dominant. 

For a typical mesoscopic silver wire of electron density 6 x 10^^ cm^"^, the low-temperature transport parameters 
are [| Lin = 1 cm > L = 30 /im and Lc\ = 50 nm. The threshold voltage is then V > 40 mV for the breakdown of 
the diffusive regime. In a two-dimensional mesoscopic channel at density 2 x 10^^ cm^^ the threshold is appreciably 
lower, with Q L = 17 /im, Li„ = 6 /im and Lei = 1-4 /im giving V ^ 0.25 mV. This shows how readily mesoscopic 
devices, particularly low-dimensional ones, can enter the non-linear regime beyond diffusive theory. 

The present is the first of three studies covering the essential formalism for non-equilibrium noise, the action of 
Coulomb correlations in non- uniform systems |p9| and finally the kinetic description of shot noise |]30| . Throughout, 
we follow one overarching principle. It is that a consistent model of non-equilibrium fluctuations will conform to basic 
properties of the electron gas in a natural way, if and only if such a model is grounded explicitly in the theory of the 
electronic Fermi liquid. 

To require that a kinetic description of fluctuations and noise respect fundamental conservation laws in equilibrium, 
is to place a unifying constraint on its low- and high-field forms together. Our aim is to catalogue all the physical 
consequences of this assertion. A viable kinetic model will necessarily recover the fluctuation-dissipation relation 
but it must also contain the equally fundamental conserving sum rules |]3^ . 

Much of the authoritative literature on diffusive noise theory, if not all of it, prefers to make a virtue of its heavy 
dependence on drift-diffusion equivalence and the fluctuation-dissipation theorem, while remaining oblivious to every 
other basic sum-rule requirement. Such an understanding is too scanty. It is the full physics of Fermi liquids which 
governs their fluctuations, a fact which cannot simply be ignored. 

An established tool for scattering-dominated noise in degenerate conductors is the semi-classical Green-function 
approach p6| , p7| , ^3| -|35|] . We take it beyond its well- understood role in time-dependent response, by proving that 
the dynamical Green function also governs the structure of the adiabatic (steady-state) fluctuations. In turn, these 
determine the mean initial strength of the time-dependent fluctuations. 

The non-equilibrium adiabatic fluctuations are exact closed functionals of their equilibrium form. This offers the 
key to practical and ffexible calculations over a wide range of useful non-pcrturbativc collision models. (Here the 
detailed classical studies of Stanton and Wilkins show the way.) Such problems are entirely out of range for 

the Boltzmann-Langevin models, of wide currency but quite narrow practicality for conductors in strong driving flelds 
|p6| , with strong internal interactions 1 36 1 . 

In section 2 we present a wholly conventional Boltzmann description of carrier fluctuations in non-uniform metallic 
systems, down to the same mesoscopic range accessible to alternative (diffusive) models HJl^. We demonstrate the 
quantitative connection between fluctuations and power dissipation well out of equilibrium. In the weak-flcld limit, 
this connection is the canonical fluctuation-dissipation theorem (FDT). At high fields, it describes the hot-electron 
contribution to current noise. This lets us calculate, in section 3, the excess thermal spectrum, which is strongly 
suppressed in a degenerate system. In section 4 we sum up, and preview the two forthcoming works. 



II. THEORY 



The theoretical discussion is in six parts. We begin by formulating the transport problem as a direct mapping of 
the electron Fermi liquid to its non-equilibrium steady state. Next we describe the steady-state fluctuations, after 
which we discuss time dependence, then the dynamic fluctuations and their formal connection with the steady state. 
This produces a self-contained expression for the current-current fluctuation, which determines thermal noise. Last, 
we analyze the connection between fluctuations and dissipation in the non-equilibrium region. 



A. Transport Model 



The semi-classical Boltzmann transport equation for the electron distribution function fa(t) = /s(r,k, t) is 



d d eE(r,t) d 

1- Vks 

dt dr h 9k 



fait) 



X [W^a'a(l - fa') fa " W^aa' (1 " fa) fa 



(1) 



Label a = {k, s, r} denotes a point in single-particle phase space, while sub-label s indexes both the discrete sub-bands 
(or valleys) of a multi- level system and the spin state. The system is acted upon by the total fleld E(r, t). We study 
single-particle scattering, with a rate Waa' = S{r — r')^/^^/ (k, k'; r) that is local in real space, independent of the 
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driving field, and that satisfies detailed balance: Wq'q(1 — = Waa'{^ — fa^)fa'^ where /^^^ is the equihbrium 

distribution. In a system with dimensions, we make the following correspondence for the identity operator: 



{l)(r)4k'} ^ Sss'Sir - r')(27r)'^<5(k - k') 



The volume ^(r) of a local cell in real space becomes the measure for spatial integration, while its inverse defines the 
scaling in wave- vector space for the local bands {k, s}. 

The first step is to construct the steady-state solution /„ = fa{t oo) explicitly from f^"^, which satisfies the 
equilibrium form of equation (|^): 



d eEo(r) d 



dr 



(2) 



The internal field Eo(r) is defined in the absence of a driving field. The quantities f'^'^ and Eq are linked self-consistently 
by the usual constitutive relations, the first being the Poisson equation 



d_ 
dr 



•eEo 



-47re 



(3a) 



in terms of the dielectric constant e(r), the electron density (/°'^(r)) = ri(r) ^X^k s/q'^' positive background 

density n+(r), which is taken to be independent of the driving field ^t\ . Normalization to the total particle number 
is J2r ^('")(/°'^('")) = The second relation is the form of the equilibrium function itself, 



J?cq _ 
J a 



1 + exp 



(3b) 



at temperature T. The conduction-band energy — es(k; r) may have structural parameters that depend on position 
implicitly. The locally defined Fermi level </'q = /i — Vq (r) is the difference of the global chemical potential /i and the 
electrostatic potential Vo(r), whose gradient is eEo(r). 

Define the difference function ga = fa — fa'^- From each side of equation (|l|) in the steady state, subtract its 
equilibrium counterpart [p8|. We obtain 



Vks 



dg^ eE(r) 9g„ e(E - Eq) df^^ 



dr 



go 



dk h 9k 

+ Y^^Wa'a - Waa'KfJ^ga + ga'f, 



Waa'ga') 

"""^gcga'). 



(4) 



The solutions to equations (^ and (Q) are determined by the asymptotic conditions in the source and drain reservoirs, 
be it at equilibrium or with an external electromotive force. The active region includes the carriers within source and 
drain terminals out to several screening lengths. This means that local fields are negligible at the interfaces with the 
reservoirs; in practice, one shorts out the fields so that E(r) — Eo(r) = beyond these boundaries. Then Gauss's 
theorem implies that the system remains globally neutral: 



^l](r)(5(r)}EE^5„ = 0. 



(5) 



Recast equation (Q) as an integro-differential equation. The inhomogcneous term on its right-hand side generates 
the explicit dependence on the equilibrium state of the system: 



(6) 



The net non-equilibrium field E — Eq is represented here as E = Eoxt + Ejnd, consisting of Eoxt(r), the externally 
applied field plus the local induced response Eind(r)- The linearized Boltzmann operator B[W^f] is 



B[W^f]. 



ak' 



Wa 



(7) 
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with W^^, = Waa' — Wa'a- Note that = if the scattering is elastic or if a hnear approximation (such as Drude) 
replaces the explicit Boltzmann collision term. 

If it is to represent the physical solution, g must vanish with E in the equilibrium limit. This is guaranteed by the 
Poisson equation for the induced field, 

^■eEind = -47re((/(r)} - {rir))) = -4^e(ff(r)). 



B. Steady-State Response 



To calculate the adiabatic response of the system about its steady non-equilibrium operating point, we introduce 
the propagator |Eg| 



dcf 5ga 



with a global constraint following directly from equation (^: 

Gaa' = for all a'. 



(8) 



(9) 



The equation for G is derived by taking variations on both sides of equation 



^BiW^fUGf,, 



eE(r') 



_d_ 
"dk' 







(10) 



The variation is restricted by excluding the reaction of the local fields Eo(r) and E(r). This means that G is a 
response function free of Coulomb screening. Here we treat the electrons as an effectively neutral Fermi liquid. In 
our second paper we will describe the complete fluctuation structure, with Coulomb effects p9[ |. 

All of the steady-state fluctuation properties induced by the thermal background will be specified in terms of G 
and the equilibrium two-body fluctuation. This consists of the "proper" electron-hole pair correlation in its static 
long- wavelength limit (up to a normalization factor). In the free-electron approximation that correlation is 



lim 



lim 



/°q(r,k-q/2) 



fcq 



(r,k+q/2) 



hcj ~ £s(k-hq/2; r) + e,(k-q/2; r) 



(11a) 



where the net momentum transfer becomes negligible relative to the Fermi wave number fcp. The formal statistical- 
mechanical deflnition of thc_(mean square) occupation-number fluctuation A/^'^ takes it as the variation of the 
occupancy of equation (3b) with respect to the electro-chemical potential and normalized to the thermal energy, 
keeping T and the local volume il{r) fixed. That is. 



A/: 



oq 



eq 



(lib) 



It is easy to derive the free-electron form of the equilibrium fluctuation: 

A/S^ = /^'^(l-/^^). 



(11c) 



When there are strong exchange-correlation interactions the two-body correlation, equation ( 11a ), is renormalized by 
a coefficient that depends on the Landau quasi-particle parameters. This coefficient carries through in Af'^'^. In the 
present model we neglect exchange-correlation effects. This is a valid approximation in dense degenerate systems 

Before discussing the non-equilibrium fluctuations we comment on the crucial contrast between the quantum- 
Fermi-liquid origin of equation ( |ll| ) and the widespread Boltzmann-Langevin approach, which is essentially classical. 
Equation (11a) describes an elementary and kinematically coupled electron- hole excitation [a fact that is particularly 
obvious in the particle-hole structure of the energy denominator, — £s(k-|-q/2; r) -I- es(k— q/2; r)]. Its form is 
determined by the same quantum dynamical equation whose semi-classical limit is the Boltzmann equation itself . 
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The physical character of a polarized excitation demands its representation as a self-contained entity. In thermo- 
dynamic terms, stochasticity attaches to the spontaneous generation of electron-hole pairs in the system and not to 
their electron and hole constituents independently (as if pairwise charge balance were of no real consequence). By 
its nature, a pair fluctuation cannot be decomposed ad hoc into two stochastically unlinked single-particle factors. 
The Boltzmann-Langevin approach, on the other hand, is tantamount to such a notional decomposition plf . Neither 
natural nor necessary for the right description of the elementary fluctuations in a charged Fermi liquid, it is eminently 
dispensable. 

Define the two-point particle-hole function A/^^l = {Saa' + Gaa')^f^- The steady-state distribution of the local 
number fluctuation is the sum of all of the two-body terms: 

A/„ = ^A/('j, ^A/^'i + J^G^^-A/^? foraUa. (12) 

a' a' 

Once the explicit solution for G is obtained, the behaviour of A/ is known completely. This non-equilibrium fluctuation 
satisfies the linearized steady-state Boltzmann equation: 

Y,B[W^fUMp=0- (13) 

/3 

We see that the solution to this equation is manifestly a linear functional of its equilibrium counterpart. Jointly, 
equations (|l2|) and ( |l3| ) mean that any scaling behaviour exhibited by the fluctuations at equilibrium must also be 
exhibited by the actual fluctuations for the non-equilibrium problem. That is the direct result of local equilibrium in 
the asymptotic leads, and of overall neutrality in the system made up of conductor plus leads. 

The proportionality of all thermally induced noise to ambient temperature T is inevitable in the degenerate limit. 
This has implications for understanding shot noise in metallic conductors |^,^,^. Shot noise never scales with T. 
We comment further on the scaling issue in the later sections. 

One may compare the analysis in terms of G and A/'^^^ with the equal-time correlator introduced by Gantsevich 
et al. [Q. The equal-time correlator is itself a hierarchical functional of other correlators (such as the current 
fluctuations), whose solutions are unknown a priori and which must be closed by force, either by truncation or by 
an ad hoc heuristic device such as Boltzmann-Langevin. This makes for a less-than-tractable computational scheme, 
at least beyond a narrow repertoire of special limits (linear response; classical particles; weak non- uniformity ) . In 
particular, a calculable strategy for degenerate non-equilibrium fluctuations, based on the equal-time correlator, has 
yet to be demonstrated. 

By contrast, we show below that G is explicitly determined by the dynamical Green function for the linearized 
Boltzmann equation [see equation (^o|)]. As pointed out by Stanton the Boltzmann-Green functions are much 
more straighforward to compute for a wide range of collision models. This ease of calculation extends to G and hence 
to A/^^\ which provides the initial conditions for a naturally closed dynamical solution. 

Global charge neutrality requires that the total fluctuation strength over the sample, AiV = ^j.f2(r)(A/(r)), be 
conserved. This constrains not only the steady-state but also the time-dependent fluctuations. 



C. Time Dependence 

Calculation of the dynamic response requires the time-dependent Green function [^3| 

R.At'tY^e{t~t')^l^ (14) 

with initial value Raa' (0) = 5aa' ■ As with G, the variation is restricted. The linearized Boltzmann equation satisfied 
by R{t — t') is derived from equation (|l]) and takes the form 



|<5„0^ + B[W^f]^p^Rfic,'{t - t') = 5{t - i')<5aa' 



(15) 



Summation over a on both sides of this equation gives zero contribution from BlW^ J]af)Ri3o.' {t - t'). Sub- 

sequent integration over t leads to conservation of normalization : 

Y,Ro,o.'{t-t') = e{t-t'). (16) 
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The time-dependent propagator is a two-point correlation. It tracks the history of a fluctuation of unit strength, 
created in state ol at time t' . The fluctuation strength in state a, at the later time t, is Raa'it — t')- In the long-time 



limit equation 
with equation 



|) goes to the steady-state equation ( [l^ ) independently of a', so that Ro 
^) this gives |Q 



Rota' {t 



oo 



AiV' 



'(t — > oo) cx A fa- Together 



(17) 



All of the time-dependent fluctuation properties induced by the thermal background are specified in terms of 
R and the steady-state non-equilibrium fluctuation A/. From the dynamical particle-hole propagator p3], that is 
^faa'i^) = Raa'{t)Afa' , One coustructs the lowest-order moment 



A/„(t) = ^A/f,W 



(18) 



in analogy with equation (^2|). Equation (|l5|), with its adjoint |Q, implies that Afa{t) = A/„ for i > 0. Thus the 
intrinsic time dependence of A/(^^(i) is not revealed through this quantity [Q. Equation ( p^ ) ensures constancy of 
the total fluctuation strength: X;r^(r)(A/(r, t)) = AiV for t > 0. 



D. Dynamic Correlations 

We move to the frequency domain. An important outcome is the quantitative link between fluctuations and resistive 
power dissipation in the non-equilibrium regime. This requires expressing both the difference function g and the 
adiabatic propagator G directly in terms of the dynamical Green function. The Fourier transform K{u!) — J dte^'^*R{t) 
of the retarded time-dependent Green function satisfies 

}R^q/(cj) = 5aa', (19) 

showing that R(a') is the resolvent for the linearized Boltzmann operator of equation (0). From equation (^6|), the 
global condition on the resolvent is 

Y,R^^,{u)^- \. (?7^0+). (20) 

At face value this does not match the corresponding criterion for G, equation (||). To solve equation ( [lO| ) for the 
steady-state propagator explicitly in terms of the dynamic one, we follow Kogan and Shul'man [p3[ in introducing the 
intrinsically correlated part of R(ct;). This is 

C..(.)=R..(.) + ^^|^. (21) 

Once the long-time adiabatic term is removed, G{lo) conveys the purely transient response of the system. It satisfies 
a pair of identities |^^. First, the Fourier transform of the relation A/(t) = 0{t)Af translates to 

Cca.{uj)AU' = for aU a, (22) 

a' 

while equation (^0|) leads to 

J2 Caa' (w) = for aU a'. (23) 

a 

The latter parallels the constraint on G. Like R(u;), the correlated propagator is analytic in the upper half- plane 
Im{w} > 0, and satisfies the Kramers-Kronig dispersion relations. Unlike R(w), however, C(ci;) is regular for w — > 0. 
We now obtain g and G in terms of the correlated dynamical response. Consider the equation 
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inversion with the resolvent yields 



eE(rO da 
h dk' 



a'f3 



(25) 



The dominant low-frequency component of R(a') does not contribute to the right-hand side of this equation. In the 
first term it results in a decoupling of the summation over a', yielding zero because i9/°?/9k' is odd in k'; in the 
second term, decoupling means that the double summation over a' and (3 vanishes by antisymmetry. In the static 
limit equation ( p^ ) becomes the inhomogeneous equation (|6|); moreover equation ( p3| ) means that g(w = 0) satisfies 
equation (m, the sum rule for g. Therefore g = g(0), or 



:c„a'(o) 



h ' dk' 



a' 13 



(26) 



This identity is central to the fluctuation-dissipation theorem. 

In models with symmetric scattering, is zero and the adiabatic Green function assumes a simple form on 
varying both sides of equation (E6|): 



n ok' 



(27) 



More generally, an analysis similar to that for g{uj) can be used directly for the adiabatic propagator. Introduce the 
operator G(ci;), defined to satisfy the dynamic extension of equation (|lO|), 



eE(r') d _^srwA 



This has the solution 



W^o,,go 



(28) 



W, 



(29) 



In the first term on the right-hand side, the low-frequency component of R(a;) makes no contribution after decoupling 
because the physical distributions on which G(u;) operates vanish sufficiently fast that '}2,^dFa/dk = 0. In the 
second right-hand term the uncorrelated parts of Kaa'ii^) and Kapico) cancel directly. We conclude as before that 



Gaa' — Cq,q'(0) ^~^'^J^ ~ E (^""'(0) ~ Cq/3(0) 



d 

dk- 



(30) 



This is a crucial result. It shows (i) that the adiabatic structure of the steady state, through G, is of one piece with the 
correlated dynamic response (the result of causality and global charge neutrality), and (ii) that the non-equilibrium 
correlation structure evolves expressly out of the equilibrium state, through the specific functional form of GAf°'^. 

We have proved that this non-perturbative kinetic description of fluctuations is self-contained, given its conventional 
set of assumptions and boundary conditions. The kinetic formalism has inherent predictive power. Hence, extraneous 
phenomenologies are not needed to make it viable. This is in sharp distinction to the diffusive Boltzmann-Langevin 
viewpoint |31|. 



E. Spectral Density 

The vehicle for the physics of current noise is the velocity auto-correlation. It is a two-point distribution in real 
space, built on the correlated part of the two-particle fluctuation Af^^, (w) — Kaa'{^)Afa' ■ Following Gantsevich, 
Gurevich, and Katilius |3^] it is a double sum over the kinematic states: 

((vv' Af (r, r'; c.)}^ 1^^^ J-^ ^ ^ Vk.Re{C„„- (.;)}vk.,, A/„, . (31) 
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Its physical meaning is the following. At any time, the system in steady state has a fluctuation background that is 
fed by spontaneous energy exchanges with the (equilibrium) thermal bath. The average strength of the fluctuations 
is fixed by the distribution A/. The elementary modes making up this background are long- wavelength electron-hole 
excitations; these are given by Af(2)(tj). The pair excitations are not themselves dynamically stable. Their transient 
evolution is determined by the propagator C(w) acting upon the ensemble-averaged background source, A/. Finally, 
the velocity-velocity correlation for the pair process is obtained by attaching velocity operators at the start and end 
of the electron-hole excitations, and summing over states [Q. 

This approach to auto-correlations makes straightforward, and completely standard, use of the Boltzmann jsj] and 
Fermi-liquid theories. Our particular contribution is to have given an explicit recipe for computing the steady- 
state form of A/ semi-classically, by analyzing the underlying adiabatic propagator G. Practical calculations should 
thereby become easier for degenerate systems at high driving fields. 

The one-point object derived from equation (pT|), 

5/(r,c.) = e^^l7(r')(((E(r).v)(E(r')-v')Af(2)(r,r»))l, (32) 

measures the local effect of fluctuations that are spread throughout the system. Formally it is the auto-correlation 
function of the power transferred from field to carriers, an inherently volume-distributed property that is represented 
here in terms of a locally defined spectral density. Sf is closely related to the thermally induced current noise, 
integrated over the entire structure. In the weak-field limit it satisfies the FDT. 
The two-point velocity correlator 

((vv'Af(2)(r,r';u;)))l/A7V, 

which is the response to a unit change of total particle number (and which does not scale with T) , should provide the 
direct basis for shot-noise calculations across distances |r — r'| comparable to the mean free path. It is natural to ask 
how shot noise fits into the framework of equation (^2|). Within semi-classical kinetics, the short answer is that shot 
noise cannot be encompassed by the generic spectrum for thermal noise. For, as we have rigorously shown, all thermal 
fluctuations are required to scale with T in the strongly degenerate (metallic) regime. Shot noise, on the other hand, 
has no such scaling. Therefore, whatever the kinetic description of shot noise may be, it is impossible for it to exhibit 
the smooth physical "cross-over" into thermal noise that is the primary feature of every diffusive model ||3l|] . 

Our approach to the kinetics of mesoscopic shot noise is explored in reference In reference we propose a 
quite specific experimental test of our theory. The new predictions made there are in stark contradistinction to the 
diffusive ones. 



F. Fluctuation and Dissipation 



The fluctuation-dissipation relation near equilibrium ties the spectral density of the thermal current fluctuations to 
the dissipative effects of the steady current in the system. However, dissipation by itself does not exhaust the physics 
of this sum rule. There are non-linear terms, negligible in linear response, that dominate the high-field behaviour of 
the noise |^,^. In view of this, it is imperative to reveal the precise nature and action of these terms. We do so. 

The resolvent property of R(ti^) provides a formal link between the steady-state (one-body) solution g and the 
dynamical (two-body) fluctuation Af^^^ at the semi-classical level. Taken to its equilibrium limit this becomes the 
familiar theorem. The connection is made in two steps. Consider the kinematic identity 

and apply it to the leading term on the right-hand side of equation (|26|). The result is 



9a 



J2Co.a'{0)iE.^)a'Arj + h^, (34) 



in which ha — J2a' i3^<^oi' {0)ga'W^, ^gp. Evaluation of the current density according to J(r) = —e{\'g), means that 
the power density P(r) = E(r)'J(r) for Joule heating can be written as 

P{r) - ^^^^(E.v)„C„„,(0)(E.v)„,A/:?-e(E.vft). (35) 
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In the second step we take the one-point spectral function Sf in the static limit, substituting for A/ from equation 
( |l^ ) in the right-hand side of equation ( ^2|) to give 



Sf {r,0) = ^^^5](E.v)„C„„,(0)(E.v)„,A/:? + 5,(r,0), (36) 

where S'g(r, 0) is generated by replacing A/ with Ag = Af — A/'^i in equation and subsequently in equation 
(^). Direct comparison of equations ( p5[ ) and ( |3^ ) leads to 

/cbJ KbJ 

This is the precise connection between the non-equilibrium thermal current fluctuations and resistive dissipation in 
the system. 

The limiting weak-field form of equation ( ^7| ) is easily obtained. We prove that it is the linear fluctuation-dissipation 
theorem. Observe that the term in h on the right-hand side varies as Eg^, while the final term varies as E^Ag; both 
of these contributions are therefore of order E^. Suppose that the system is uniform. Then E = Eoxt = E acts along 
the X-axis. Division by E'^ on both sides of equation ( ^7|) gives 

(38) 



E^ knT E 



where a is the low-field conductivity. Equation (|3^) is the canonical FDT. 

The non-dissipative and purely non-equilibrium structures beyond P(r) can be expanded similarly to it. We discuss 
the symmetric-scattering case, for which there is no contribution e(E'v/i). Within Sg we apply the formula for the 
adiabatic propagator, equation (p7|), to express A;; ~ ^ GA f°'^ in terms of the correlated dynamic response function 
C(a;). This produces the closed form 



^.('■'O) = j4EE(E-v)oC.,(0)(E.v),(^C,.,(0)^-^' 



J- 5] 5:(E.v)„(C(0)E.v)L,(1 - 2rj)AfJ. (39) 



kBTQir) ^ , 

k,s a' 

The second line follows from the first after using equation ( |3^ ) to express dAf°'^/dh in terms of and Af'^'^, and 
taking an inner sum into (C(O)E'v)^. 

The expression above differs markedly from the rate of energy loss P{r) by Joule heating. In contrast, Sg{r,0) 
relates directly to non-equilibrium broadening of the fluctuations, due to the excess energy gained from the field 
during intervals of ballistic flight ||2^ ] . The extent of the broadening is limited by dynamical dissipation of the excess 
energy, locally (by prompt inelastic scattering) or remotely (by carrier relaxation in the ideally absorbing terminals) . 
The impact of this term on current noise is felt only for substantial departures from the weak-field regime. 

There exist several alternative generalizations of the FDT for extended bulk systems ]4^-|4^. We mention the 
best known, which defines the non-equilibrium noise temperature Tn pivotal to the interpretation of device-noise data 
p^ . Phenomenologically Tn is obtained, for a non-linear operating point, by normalizing Sf with the differential 
conductivity (Jx{E) = dJx/dEx such that kBTn{E) = Sf/(Tx{E)El, corresponding to the output of a small-signal 



noise measurement. (In general Tn is not isotropic.) Our equations ( |35D - (39) provide a microscopic framework for 
computing the noise spectral density in a wide class of degenerate systems. Since <Tx{E) is also calculable within the 
same framework, this yields T„. 



III. APPLICATION TO HIGH-FIELD NOISE 



We can now explore one of the most significant properties of the excess spectrum Sg: its strong inhibition by 
degeneracy. That there exists an additional, purely quantum-statistical, constraint on field-driven broadening is seen 
directly in the factor (1 — 2/°'^) of equation (p9|). This suppresses the contribution of Sg relative to the corresponding 
classical result, in which the factor is unity. Suppression of electron heating by Pauli exclusion refiects the large energy 
cost of displacing electrons deep inside the Fermi sea. 
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To highlight the difference between dissipative and hot-electron terms, we revisit a simple example |3^,Q, the 
uniform electron gas in the constant-collision-time (Drude) approximation subject to a field E = —Ex. Expressions 
for the power density P and hot-electron component Sg are derived in the Appendix. The thermallydriven current- 
current spectral density, taken over a uniform sample of length and total volume fl, is given by ||44|] 



siE, 4 y: m E (-5^) (-^) Af(^) (-) 



Writing the sample conductance as Q = QP / L^E^ , the static limit of the spectrum is determined by equation (|37|): 



(40) 



1 



^^(O) 



PknT 



= ^GkBT 



1 



An / m*filE^ 



(41) 



We have substituted for P and Sg respectively from equations (A7) and (A9). The electronic density is n while 
An = AN/fl is the number-fluctuation density. The effective electron mass is m* and fi^ is the mobility. 

The term Sg/PkBT is a relative measure of the hot-electron contribution to the noise. The inhibiting effect of 
degeneracy, through An/n, is greatest at low temperature and least in the classical regime. When the Fermi energy 
Ef satisfies ep ^ kBT, the ratio An/n goes to unity and the hot-electron term is that of a classical electron gas (low 
density, high temperature). Its form in the high- field limit E y^ksT/m*^'^ is 5 ^ 4:Qm*^lE^, asymptotically 
independent of T. 

On the other hand, when kBT <C ep the system is strongly degenerate. In a i/-dimensional system we have ep oc n^^"^. 
Then 

An fceT" dn vkBT 



with equation (^) this leads to 



n n dep 2e 



SiE,0) ,^^./nO^, ^^^^ 



5(0,0) 2 V £f 



Note that the thermal fluctuation spectrum S{E, 0) necessarily vanishes with temperature, while its ratio with the 
Johnson-Nyquist spectral density 5(0, 0) — AQkBT continues to exhibit a hot-electron excess which is now scaled by 
the Fermi energy. 

Figure 1 displays the excess-noise spectral ratio in a two-dimensional electron gas, as a function of the applied field, 
when T ranges from the degenerate limit to well above the Fermi temperature Tp = ep/fce. For T much greater 
than both Tp and m* ji^E'^ / kB the excess contribution becomes classical, independent of temperature, and thus small 
compared with the now-dominant base value 5(0, 0). This is evident in figure 1 through the gradual downward shift 
of the plots, with increasing T . 

Equation ( ^ ) may be compared with a perturbative estimate by Landauer |47[ in the degenerate limit, for which 
the analogous excess term is {5U/kBT)^, where 5U ~ m* fjLcEvp is a characteristic energy gain. Taken at face value, 
this would suggest that hot-electron effects in the low-T regime can be enhanced even more by further cooling of the 
system. 

This counter-intuitive result comes from inappropriate use of perturbation analysis. Series expansion of the thermal 
current noise, in powers of E, fails to account for non-analyticity of the full non-perturbative solution in its approach to 
equilibrium ||48|] . Non-analyticity of the distribution function /k precludes the reliable calculation of moment averages 
by expanding about equilibrium, as in reference |47| . (It is reassuring - and only seemingly fortuitous - that the 
actual linear current response is reproduced exactly by solving the transport equation, as usual, to first order in the 

field m.) . n 

The relevance of non-analyticity to transport physics has been questioned by Kubo, Toda, and Hashitsume [[49| . 
They regard its appearance as spurious, a specific artefact of the crude way in which the Drude approximation treats 
real collisions. That is to overlook the appreciably broader evidence for non-analyticity in the variety of collision 
models assembled by Bakshi and Gross p^ . 

Even within the Drude model of a degenerate conductor (over-simplified though it is), it is clear that its exact non- 
perturbative solution does produce physically consistent scaling of the excess noise with T. Equally clearly, finite-order 
response theory does not. Kubo linear response recovers only 5(0, 0) and misses the non-linear excess noise altogether. 
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Such sharp differences between perturbative and non-perturbative predictions should be experimentally measurable 
in the hot-electron spectrum. In our view, issues of non-analyticity and its physical manifestation remain open. 

We make some final comments on shot noise and the impossibility |2^j2l|,^ of a theoretical cross-over, unifying 
thermal and shot noise for mesoscopic metallic wires. The diffusive cross-over formula, ostensibly identical in form 
to equation (EO), always generates a T-independent term p4-pi|. One might have expected that a computation of 



the spectral density of equation (40), taken in the semi-classical quasi-ballistic limit <C m*VFiic/e, would yield an 



expression for S{E,0) that is independent of T and proportional to the current / = QV; in other words, shot noise 



We have carried out this quasi-ballistic exercise for a degenerate system, in simplified form |4^]. At high fields it 
gives S{E,0) ~ 2eI{k^T / e^). This is indeed linear in / but thermal nevertheless, since its immediate source is the 
generic spectral relation, equation (|39|). 

Thermal fluctuations are induced by spontaneous and quasi-continuous changes in the total internal energy of 
carriers, throughout the whole active volume of a device. Shot-noise fluctuations are induced by spontaneous and 
discrete changes in total carrier number, through the device's interfaces with the outer circuit. Such qualitative and 
topological distinctions may be of little practical importance in the classical macroscopic world. However, it is not at 
all clear that they are immaterial to the metallic mesoscopic regime. The issue is under active examination 

A leading task is to identify the kinetic origin of the empirical cross-over between thermal and shot noise, apparent 
in real mesoscopic conductors Once again we stress that, regardless of how shot noise is to be described 

microscopically, the logical and conceptual gaps between diffusive explanations of the cross-over (quantum as well 
as semi-classical) P, p7| , ^ and strictly conventional kinetic theory have already been uncovered, characterized and 
analyzed |20,El[. We will present a fully detailed semi-classical kinetic model of shot noise in due course. 



IV. SUMMARY 



We have described, and applied, a genuinely non-equilibrium kinetic formalism for current fluctuations. It holds for 
metallic systems down to mesoscopic scales, within the ambit of semi-classical theory. Our strategy for incorporating 
microscopic Fermi-liquid correlations within the Boltzmann picture safeguards the conservation laws at the two-body 
level. Conservation continues to underpin the nature of current noise at high flclds. 

Our theory leads to a precise quantitative link between non-equilibrium thermal current fluctuations and energy 
dissipation. In its low-field form, this is the standard fluctuation-dissipation theorem of linear-response analysis. At 
high fields, it highlights the pervasiveness of strong degeneracy even in hot-electron noise. 

In a completely standard description, such as ours, correctness and calculability do not issue from the ad hoc 
assumption of fictive Langevin noise sources; even less do they rely on diffusive analogies that fail manifestly to respect 
the canonical sum rules in electronic Fermi systems. Rather, the model's integrity will stem from the microscopic 
structure of its underlying Green functions. They describe how the fundamental electron-hole pair excitations evolve 
within a metallic conductor, in the semi-classical limit. 

We have discussed how to map these native polarized correlations non-perturbatively, from their equilibrium dis- 
tribution to its analogue in the externally driven conductor. The resulting high-fleld noise spectrum yields a faithful 
signature of its source: the elementary non-equilibrium electron-hole polarization processes. 

The main, and physically inevitable, consequence of this rigidly orthodox kinetic investigation is the intrinsic 
scaling of degenerate-electron fluctuations with thermodynamic temperature. For diffusive phenomenologies, this is 
one phenomenon too many. To sustain their predictions for the shot noise of metallic wires, they have no choice but 
to deny outright all possibility of T-scaling for hot-electron noise |^ . This must be so for any theory that predicts a 
seamless cross-over between thermal and shot noise. 

There is an undeniable connection between T-scaling and the dominance of long-range screening, as of degeneracy, 
in the polarizable electron gas. Logical examination shows that this nexus can be broken only by contradicting the 
standard picture of charge fluctuations in metals. Diffusively inspired models would seem to do exactly that [p^|-pl|. 
So far, no such model has rationalized the heroic departure from principles that have been understood, widely and 
thoroughly, for some time ||3^ ]. 

We envisage two extensions to this work: the systematic inclusion of Coulomb screening within the microscopic 
structure of the fluctuations p9| , and the analysis of shot noise as a kinetic process quite separate from thermally 
driven noise pO[ . Coulomb effects are particularly evident in strongly conflned electron systems, such as the two- 
dimensional electron gas in a III-V heterojunction quantum well |2^ . Self-consistent Coulomb screening in a conflned 
channel should markedly reduce the scale of thermal fluctuations in the current. 

Shot noise and thermal noise have disparate properties, which no-one disputes. Shot noise never scales with ambient 
temperature, while excess thermal current noise must do so if there is strong degeneracy. Coulomb effects too may 
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differentiate between the two kinds of fluctuations. If so, then selective action of the Coulomb correlations could 
serve as an experimental tool to distinguish between excess thermal current noise and shot noise. This would help to 
pin-point both the distinct sources of non-equilibrium mesoscopic fluctuations and the disposition of Coulomb forces 
at small scales. We take up these themes in the forthcoming papers. 

Every formalism for mesoscopic noise stands or falls by its new predictions. Ours is no exception ||2^,|2^. Boltzman- 
nian kinetics are obviously not equipped to give the final word on quantum fluctuation effects; be that as it may, it 
hardly needs saying that any mesoscopic model, whatever its origin, should be totally consistent with the established 
physical facts. In the context of the metallic electron gas, noise descriptions which claim to be truly microscopic must 
address full sum-rule consistency as a matter of course. This transcends semi-classical analysis and is by far our most 
important message. 
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APPENDIX A: UNIFORM DRUDE MODEL 



We derive the dynamical fluctuation structure for a single parabolic conduction band with uniform electron density 
n and constant mobility fi^ — eT/m*, where r is the spin-independent collision time and m* the effective mass. The 
system is driven by a uniform field E = E — i?x acting in the negative (drain to source) direction. We take variations 
which are homogeneous over the sample region, so that the fluctuations of interest have no spatial dependence. 

The Boltzmann equation in the model is 



d eE d 
dt h dkx 



hit) = 



(fit)) fT 
in) r 



(Al) 



Since the Boltzmann operator is linear, the fluctuation structure is qualitatively similar to that for elastic scattering 
[differences arise from the inhomogeneous term in /°'^, notably in the behaviours of R{t) and A/(t)]. We solve equation 
(Al) by Fourier transforms in reciprocal space, so that the transform Fp = ri^-'^^^./k exp(ik-p) of the steady-state 
distribution takes the form 



^ ~ FS"" 1 



ikdPx ' 



(A2) 



where fc^; = eEr/h and Fq — \{f) per spin state. While a formal distinction is made between Fq and Fq'^, the 
physical normalization is always Fq — F^'^ = ^n. Note also that Fp is singular for p^. = —ik~^^. In wave-vector space 
this means that /k is non-analytic at i? = 0. The same can be said for A/k. 
The transform of the dynamic response function, 



■Rpp, [uj) 



1 

f72 



has the equation 



EkEk'I^kk'(c.') exp[z(k.p - k'.p')], 



[-iujT - ikdPx + 1] TZpp'iuj) ^ Td{p - p') 



F^^ 



P • 



(A3) 



For p = this leads to 



7^op' (w) 



Sip') 



+ irj) 



(A4) 



On the other hand, the low-frequency adiabatic part of TZppi scales with the steady-state solution Fp [in a collision- 
time model the asymptotic form Fp/^n replaces AF p/ ^An]. On denoting the correlated part by Cpp' and recalling 
that the adiabatic part exhausts the normalization of TZop' , we obtain 



12 



i{uj + irj) Fq 



(A5) 



When the above is put together with equations (A2)-(A4) we arrive, after some algebra, at the expUcit formula for 
the correlated propagator: 



1 - ikdPx ~ i'^T 



(A6) 



We can use equation ( |A6D directly to evaluate both dissipative and non-dissipative contributions to the noise. Using 
the reciprocal-space representation v ^ —i{h/m*)d/dp^ the power density P of equation (BSh is 



P = 2- 



k^T \ m* 



ksT \m 



dpi " 



(A7) 



The Drude conductivity a = nepo appears when we apply the relation 

a2 



to the middle line of the equation. A contribution containing {vr^Af^^) = vanishes trivially. 

The hot-electron spectral density Sg in the static limit [recall equation (|3^)] is calculated similarly: 



e'E^T^h 



<fp' j <Fp"vxCpp, {QKCp,p„ m-tp'^AFp) 

dpa: 1 - ikdpx \dpx 1 - ikdPx )\p^Q 



" d Fp/Fo 




" d 


-¥x^Fp 


dpx 1 - ikdPx_ 


p^O 


dp'x 


1 - ikdp'x 



(A8) 



We evaluate this with the help of the relations AFq'^ — ^An and {dFp/dpx}p^o — ikdFo, the latter following from 
equation (|A2|). The result is 



,g - am p^E [ — 



At 



(A9) 
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FIG. 1. Zero-frequency spectral density of the excess (hot-electron) thermal noise, above the equilibrium noise, in a degen- 
erate and uniform two-dimensional electron gas. The excess is plotted as its ratio with the zero-field noise S{E=0) = AQk^T as 
a function of driving field E and for temperatures T between and 900 K, in steps of 150 K. The dot-dashed lino is at T = 300 
K. In the degenerate limit T <C Tf, thermal noise scales with T; thus the excess-noise ratio is independent of temperature. At 
high temperature, the excess ratio for a given field value diminishes as its denominator <S(0) becomes dominant. 
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